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Abstract 

We estimate the isospin-violating QED radiative corrections to the charged-to- 
neutral ratios of the decay rates for and in non-leptonic B meson decays. In 
particular, these corrections are potentially important for precision measurement of 
the charged-to-neutral production ratio of B meson in e~^e~ annihilation. We calcu- 
late explicitly the QED corrections to the ratios of two different types of decay rates 
T{B+ J/i;K+)/T{B° J/ipK^) and T{B+ D^D^)/r{B^ taking 
into account the form factors of the mesons based on the vector meson dominance 
model, and compare them with the results obtained for the point-like mesons. 



At present 5-factory experiments, a measurement of the B^B~ to B^B^ production ratio 
in e^e~ annihilation is essential for accurate determination of standard model parameters 
such as quark- mixing matrix elements. The B mesons are produced in pairs at the T(4S') 
resonance, and the charged-to-neutral production ratio given by 

V (T(4S) B^B^) ^ ' 

has been recently measured by CLEO 1], BABAR |2^, and Belle The experimental 
values of R^^^ typically range from 1.01 to 1.10 and the latest significant measurement gives 
= 1.006 ± 0.036 ± 0.031 ^} The theoretical prediction of the deviation 5R+I^ has 
been first made by calculating the Coulomb interaction between point-like mesons [Sj, and 
reexamined later by considering the effects of the electromagnetic form factors of the mesons 
and their interaction vertex jH]. The theoretical values of the ratio W'f^ estimated so far 
range from 1.03 to 1.25. As well as the Coulomb interactions, it was also shown that the 
strong interaction phase in the region of a strong resonance can modify the experimental 
result of the ratio near threshold in e^e~ annihilation 

The practical measurements of the ratio use the yields of events for the exclusive 
processes where the ratio of the decay rates for B^ and 

is determined by the isotopic invariance. An example of such isospin-related decays is pro- 
vided by 5+ J/i)K+ and 5° J/ipK^ where the deviation S+/^{J/ipK) due to the 
isospin-violating mass differences between these two decays is about 0.0014. Since this 
isospin violation is very small and doesn't exceed the present experimental accuracy, it has 
been assumed that r{B^ J/ipK'^) = r(i?° J/ipK^). In decays of this type, however, 
there is still an isospin violation due to QED radiative corrections. At a non-zero deviation 
5+/", the actual value of the production ratio should be corrected to _R+/° = Rtlpr/ (1 + 5"^^°). 
For precision measurement in present and future experiments, it is important to understand 
the production ratio of specific exclusive modes with accuracy better than 1%. Espe- 
cially, the decay mode such as B^ J/ ipK^ requires a better understanding due to its great 
importance in the measurement of CP violation. Hence, the deviation 5^^'^{X) for a specific 
decay mode needs to be evaluated more carefully. 

^It can also be mentioned that BABAR recently reported their preliminary result of the direct measure- 
ment of the branching fraction S(e"'"e^ — > B'^B^) = 0.486 ± 0.010 ± 0.009 using partial reconstruction of the 
decay B^ D*+l-fii g). 
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The isospin violation in B ^ D semileptonic decays was discussed earlier in Ref. jH]. In 

the decays B~ D^W and D'^W, for example, the kinematic difference due to the 

isotopic mass splitting of the D mesons in the spectra is not negligible and can amount to 
a sizeable fraction of 1%, while the difference is compensated in the total decay rate with 
a much better accuracy as predicted by the heavy quark effective theory (HQET). On the 
other hand, in the nonleptonic decays where isospin violation due to the mass difference 
between the charged and neutral mesons is expected to be very small, the QED corrections 
may not be negligible and dominantly lead to the ratio of the corresponding decay rates being 
off unity. In the standard model, the QED corrections usually amount to an a order where 
a is the fine-structure constant and differ by sizeable amounts depending on the final state 
masses. In this paper, we mainly concentrate on the two different types of the non-leptonic 
B decays, B+ J/iIjK+ and 5° J/ipK^ {VP type) decays as well as B+ D+D^ 
and B^ D^D^ {PP type) decays where PiV) denotes a pseudoscalar(vector) meson, and 
evaluate explicitly the QED radiative corrections to the ratios of those decay rates. 

J/l|r J/i)/ JA|/ J/\|; J/I|f 
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Figure 1: One- loop diagrams for B^ —>■ J/ipK^ decays. 




Figure 2: One- loop diagrams for B^ D^D^ and B^ decays. The self-energy 

makes no contribution to the ratio of the corresponding decay rates. 

The calculation of the QED radiative corrections due to virtual- and real-photon in the 
leading order for point-like charged scalar particles is quite straightforward. The radiative 
corrections due to virtual-photons can be obtained by performing the loop integrals shown 
in Fig. n and Fig. |21 However, one cannot simply treat the mesons as point-like particles 
because the QED corrections are sensitive to the structure of the mesons, and thus we require 
models for the form factors of mesons. The electromagnetic form factor Fx of a meson X 
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is successfully modeled in the intermediate energy region using the vector meson dominance 
model (VMDM) jH]. The most commonly used VMDM assumes that all photon-hadron 
coupling is mediated by vector mesons and the form factors are expressed through their 
propagators. Taking account of the electric charges of the quarks comprising the mesons, we 
use the VMDM based parametrization for the form factors in such a way that each meson 
behaves like a point-like particle in a zero photon momentum limit as follows:^ 



3 3 — k'^ 

I I ml 

3 3 — 

FKoi^) ^ FBOKOj/^k""), (3) 

2 1 ml 



3 3 — ' 
^ /,2n 2 2 ml 
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where k is the momentum of photon, and mp is the mass of p meson. Then one can obtain 
the leading virtual-photon corrections to the ratios of the decay rates by replacing the photon 
propagator l/k"^ with the modified propagators as follows: 



L 



k'^ //i R, „ \k'^ k'^ — mi 3 dmik'^ 



k'^ R„ \k'^ k"^ 



Bz ^ JBz ^ - "V/ 

A f ^,^ (1 1 ral a 



k"^ J yk"^ k"^ — m^ 3 dm^ k"^ — 

where the subscripts Ai and Bj denote each loop shown in Fig. ^and Fig. |21 Note that it is 

clear from the above equations that UV divergences arising from the one-loop calculation are 

naturally eliminated. For actual calculation, it is convenient to introduce a small non-zero 

photon mass m-^ and obtain the analytic expressions of the virtual-photon corrections for 

^We assume here that the form factors are primarily associated with the p and oj mesons, and neglect the 
(f> meson contributions for simplicity, which does not change our result significantly. In other words, we take 
into account the form factors for the electromagnetic current of the u and d quarks (2/3)u7^u — {l/3)d'j^d, 
and treat the current of the s, c, and b quarks as having a point-like form factor. 
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the point-like mesons first, so that the additional terms due to the p meson propagator in 
the VMDM can be simply obtained by replacing the photon mass with the p meson mass. 
Then, the IR divergences arising from the one-loop diagrams are expressed by Inm^, and 
the linear terms of can be dropped whereas those of irip cannot because the p meson is 
heavier than the K meson. 

After a straightforward but somewhat lengthy calculation, the leading virtual-photon 
corrections to the ratios of the decay rates are given by 



where 



fK{x,y;w) 
gK{x,y,v;w) 



§^ {fK{rj/i„ Tk; r^) + 0(rx)) 
i^{9K{rj/^,rK,rp;r^) + 0{rp)) 



(point-like mesons) 
(VMDM) 

(point-like mesons) 
(VMDM) 



-4 + 2(2 



fDix;w) 



gD[x,v;w) 




(5) 
(6) 



(7) 



(8) 



(9) 



(10) 



and where rp = mp/m^ and we use mp+ = mpo = mp. In the above equations, Inrx term 
in Eq. ((Tj) has a mass singularity in the limit rrix — > and it is removed by the same term 
in the real-photon corrections. In the VMDM there are additional terms such as mp/rriK 
for non-zero rup and rriK in general, and Eq. (jH)) was obtained for the K meson satisfying 
mp/2 < rriK < Trip. In the limit rriK — ^ 0, the additional singular terms In mi^- and mp/rriK 
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in Eq. (jH)) produced by the p meson propagator are replaced by Inr^, and the function qk 
can be rewritten as 



3 57r2 
'4^6 
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As for the D mesons, one can see the term mp/rriD in Eq- (fTIHl and this formula is valid for 
any other scalar meson P satisfying the condition mp > rrip. 

The real-photon corrections can be simply obtained by assuming that all mesons behave 
like point-like particles. 



a 
2^ 
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where 
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Now one can obtain the total QED radiative corrections ^q^d = Kirt + Klaii ^^"^ explicitly 
see that the IR divergent terms expressed by Inr^ in the virtual-photon corrections cancel 
out the same terms arising in the real-photon corrections in a usual way. 

Using the standard values of the meson masses, we estimate the numerical values of the 
QED radiative corrections to the ratios of the decay rates: 

—0.0012 (point-like mesons) 
0.0014 (VMDM) 

0.0115 (point-like mesons) 
-0.0013 (VMDM) 



X 



(14) 



(15) 
(16) 
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As a comparison, we list the radiative corrections for the point-hke mesons too. In B 
J / ipK decays, the QED radiative correction to the ratio of the decay rates is small because 
the real-photon correction is offset by the virtual-photon correction. On the other hand, the 
QED correction in B ^ D^D decays is sizable for the point-like mesons but significantly 
reduced in the VMDM due to the similar compensation between the real- and virtual-photon 
correction. In both cases, the deviations 5+^° due to the isotopic mass splitting of the K and 
D mesons are about 0.0014, so that the total isospin violation effect in S — > J/ipK decays 
becomes about 0.3% in the VMDM while that in B ^ D^D decays becomes negligible. In 
different decay modes, of course, the corrections can be enhanced or reduced due to different 
final state masses. Therefore, in order to obtain the production ratio R^/'^ more accurately 
in precision measurements, one should include the correction 5"*"/° or find an appropriate 
decay mode where the isotopic violation due to the QED effects is very small. As well as 
the determination of R~^^^, the obtained results are also of importance in the study of the 
structure of hadron. We explicitly show in Eq. ()15|) and Eq. (fT^ how the form factors of 
the mesons affect the radiative corrections. This work can be used to test the model itself 
and will be useful for the study of the form factors when future experiments are available to 
access the ratios of the decay rates. 

The author would like to thank M.B. Voloshin for suggesting this problem and for helpful 
discussions, and S. Rudaz for his useful comments. This work is supported in part by the 
DOE grant DE-FG02-94ER40823. 
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